Motion of vortex lines in nonlinear wave mechanics 
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We extend our previous analysis of the motion of vortex lines [I. Bialynicki-Birula, Z. Bialynicka- 
Birula and C. Sliwa, Phys. Rev. A 61, 032110 (2000)] from linear to a nonlinear Schrodinger 
equation with harmonic forces. We also argue that under certain conditions the influence of the 
contact nonlinearity on the motion of vortex lines is negligible. The present analysis adds new 
weight to our previous conjecture that the topological features of vortex dynamics are to a large 
extent universal. 



I. INTRODUCTION 



II. NONLINEAR SCHRODINGER EQUATION 
WITH HARMONIC FORCES 



Owing to recent advances in the experimental studies 
of Bose-Einstein condensation (BEG), a detailed theo- 
retical description of the motion of vortex lines in the 
condensate may soon become testable. In this report 
we continue our previous analysis |^] of the evolution of 
vortex lines as determined by the quantum mechanical 
wave equations. We were able to extend our methods of 
obtaining analytic solutions to a nonlinear Schrodinger 
equation. 

The standard tool in the study of BEG is the Gross- 
Pitaevskii (GP) equation that describes zero range in- 
teractions of the condensate particles. All studies of the 
condensate dynamics based on this equation must rely 
either on numerical methods or on approximations [Q, 
since analytic solutions of the GP equation in three di- 
mensions are not available. In order to obtain a rich 
family of explicit solutions allowing for elaborate vortex 
line structure, we modify the form of nonlinearity replac- 
ing the zero range interaction by harmonic forces. We 
are fully aware that atoms do not interact via harmonic 
forces but the aim of our study was to answer the ques- 
tion: Does a nonlinear term modify in an essential way 
our results obtained for the linear Schrodinger equation? 
In Section II we prove that the answer to this question 
is negative in the case of nonlinearity corresponding to 
harmonic forces. In the other extreme case of contact 
interaction leading to the GP equation we do not have 
analytic solutions to argue our case. However, the non- 
linear term in this equation always vanishes on vortex 
lines. In Section III we show that in the vicinity of the 
vortex line it remains much smaller than the kinetic en- 
ergy term. Therefore, we can argue that also in the case 
of contact interactions the effect of the nonlinearity on 
the qualitative features of the motion of vortex lines dy- 
namics at close approach is negligible. 



The starting point of our study is the nonlinear 
Schrodinger equation describing the dynamics of atoms 
in a harmonic trap interacting via harmonic forces in the 
Hartree (mean field) approximation. We shall assume 
that the trap is fully anisotropic with the characteristic 
frequencies denoted by Cbx,Cjy and Cj^- The interatomic 
harmonic forces can be either repulsive (Sl^ > 0) or at- 
tractive (il^ < 0). In natural units {h — 1, m — 1), our 
nonlinear Schrodinger equation reads 

- ^ydV|^(r',i)|2(r'-r)2V.(r,<). (1) 

Following our earlier work on systems with harmonic 
forces we introduce the following set of time de- 

pendent global quantities {N is time independent) 



jd?r\^{v,t)\\ 
RAt) = Jd^rx\^P{r,t)\'', 

PAt) = -i /dVV*(r,i)V,V(r,t), 
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and analogously for the components in the y and z di- 
rections. With the use of the quantities R(t) and U{t)^ 
Eq. (nj) can be rewritten in the form 



^{r, t) = e'fW+^'^W-'^ 0(r - b(t), t). 



(22) 



idtip{v,t) 
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^{r,t) (8) 



where = (I;2 _ ATj^a^ ^2 ^ ^2 _ ^^2^ ujI=CjI- NQ^ 
are the squared frequencies modified by the mutual in- 
teraction and U{t) — Uxx{t) + Uyy{t) + Uzz{t)- 

As a result of Eq. (|l|), the global quantities obey a 
set of ordinary linear differential equations (again we list 
only the x-components) 



dN{t)/dt = 0, 
dRx{t)/dt = Px{t) 
dPx(t)/dt = -il^lRxit), 

dUxccit)/dt = 2WxAt), 
dTxx{t)/dt = -2ujlWxx{t), 
dWxx{t)/dt = Txx{t) - ujIUxS)- 



The solutions of these equations are 
N(t) = TV, 



(9) 
(10) 

(11) 
(12) 
(13) 
(14) 



(15) 
(16) 



Px{t) = -Rx{Q)CJx sin{Coxt) + P:.(0) cos(Coxt), (17) 

, ,sin^(cL'rO . , sin(2ix'ri) 

TxxiO) ^-^ + Wxx{Q)^^. 18 



,2 ■„2t 



UJx 



Txx{t) = C/:r2:(0)w^ sin {uJxt) 

+ TxxiO) cos'^iuxt) - WxxiO)uJx sm{2u;xt), (19) 

w TT /r,N^a;Sin(2w2;t) 
Wxx{t) = -Uxx{()) ^ 

+ Txx[0 f-^^^^ + Wxx{Q) cos(2w,i). (20) 

ZUJx 

The center of mass variables R(t) and P(t) oscillate with 
the original trap frequencies {Cjx,LOy,Lbz) and the vari- 
ables U, T, and W that describe the motion of the inter- 
nal quadrupole moment oscillate with the doubled mod- 
ified frequencies {uJx,Ltjy,ujz)- 

We have found that each solution ?/;(r, t) of the nonlin- 
ear equation (|l]) is related in a simple way to the corre- 
sponding solution (/)(r, t) of the linear Schrodinger equa- 
tion 

^mr,t)=\--A + ^ f —U{r,t). (21) 

In order to prove this assertion we seek the solution of 
the nonlinear equation in the form 



Upon substituting ( p^ ) into (||) we find that (/)(r, t) must 
obey Eq.(pTl) while the scalar function f{t) and the two 
vector functions a(i) and h{t) must obey the following 
ordinary differential equations 



dsL{t)/dt = VpTl{t) - b(i), 
dh{t)/dt = a(t), 

df{t)/dt^]-{u:lhx{tf +u:lhy{r)^ 



(23) 
(24) 



&{t)-&{t) ~ VL\U{t)+'R.{t)-~R{t)/N). (25) 
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The solutions of these equations for ax{t) and bx{t) are 



ax(t) 



1 r 



-Rx{0){uOx sm{Ldxt) - LOx svaiuxt)) 



+ Px{Q){cos{ujxt) - cos{uJxt)) 



bx{t) = RxiO){cos{u}xt) - COs{LJxt)) 
+ p (o)( Sin((I'xt) _ sinjujxt) ^ 

UJx LOx 



(26) 



(27) 



and similarly for the remaining components. The phase 
f{t) can be obtained from ( p5| ) after performing a 
straightforward integration of products of trigonometric 
functions. Next, we note that the vanishing of a, b, and 
/ at i = implies that the initial values V'(r,t = 0) 
and (/)(r, t = 0) coincide. Therefore, we can argue that 
if 0(r, t) is a solution of the initial value problem for 
the linear Schrodinger equation, then 'ip{v,t) given by 
Eq. (|2^) will be the solution of the initial value prob- 
lem for the nonlinear Schrodinger equation. Note that 
the initial values of the wave function for the non- 
linear problem enter not only through the solution of 
the linear problem but also through the initial values 
R(0), P(0), C/(0), T(0), 1^(0) of aU global variables. 

The form (^2|) of the solution of the nonlinear prob- 
lem shows that there are three effects of the interatomic 
interactions. The wave function acquires a time- and 
space-dependent phase factor, the frequencies of the trap 
are modified (Cj ^ lo), and the whole wave function un- 
dergoes an additional rigid motion described by the vec- 
tor h{t). This motion is made of oscillations with the 
frequencies present in the problem. The amplitudes of 
these oscillations are determined by the initial values of 
the center of mass position R(0) and momentum P(0). 

The relation between the general solution of the 
nonlinear and the linear Schrodinger equations obtained 
in our model with harmonic interatomic forces, enables 
one to analyze the effect of this particular type of non- 
linearity on the motion of vortex lines. The phase factor 
exp(i/(t) -I- ia{t) ■ r) does not have any influence on the 
vortex lines since it never vanishes. Thus, the only ef- 
fect (apart from the modification of the trap frequencies) 
is the shift of the argument in Eq. (p2|) by h[t). This 
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time-dependent shift causes oscillations of the same vor- 
tex structure that is already present in the wave function 
4>{r, t) of the linear problem. Thus, the topological struc- 
ture of vortex lines is not affected by the interatomic, 
harmonic interactions. Vortex creations, annihilations, 
reconnections, and switchovers will occur unimpeded as 
in the linear case. 

Eq. (H^) can also serve as a practical tool to generate 
explicit solutions of the nonlinear Schrodinger equation 
(|l|) with (almost) arbitrary vortex structures. To this end 
we only need to obtain such solutions of the linear equa- 
tion. This can be done with the help of the generating 
functions, as described in Ref. 

As a simple illustration of this general analysis, let us 
consider the motion of a single straight vortex line along 
the 2;-axis displaced at t = by a from the origin in the 
x-direction. We assume that the initial wave function 
has the form 

l^vort (r, 0) 

= Nq(x - a + ly) eyii>{ 1 ), (28) 

where iVo is the normalization constant and we have cho- 
sen the ground state as the "background" wave function. 
The time dependent solution of the linear Schrodinger 
equation satisfying at t = the initial condition ( p8| ) is 
for all values of t given by 



0vort(r,O ^ iVoCXp(-«^ ^ '-t) 



(29) 



X (xe 



— a + iye 



')exp(^ 



One can obtain this solution using the general method 
introduced in Ref. . The position of the vortex line in 
the xj/-plane is determined by the zeros of the equation 



\xe ^ — a + lye "I = 0, 



whose solutions are 



X — —acos{ujyt)/ C0s((ijJ2^ — ^y)t)j 
y = -a sin(w^i)/ cos((cjj; - ujy)t). 



(30) 



(31) 



In the isotropic trap the zeros follow a circle, in all other 
cases they follow an unbounded curve; the vortex line 
moves with an unlimited speed. In order to obtain the 
solution of the nonlinear equation we need the displace- 
ment vector b(t). The only nonvanishing components of 
R(0) and P(0) in this case are: Rx{0) = a and Py{0) = a 
and the vector h{t) has only the x- and y-components 



87ra|?/'pV 



bx{t) = a{cos{ujxt) — cos{uJxt)) 



(32) 



The vector b(t) will draw a Lissajous figure in the xy- 
plane. The motion of the vortex line in the nonlinear case 
in the anisotropic trap will be the composition of an un- 
bounded motion determined by the condition ( ^ ) with 
the motion following the Lissajous figure determined by 
Eq. (il). 



III. THE ROLE OF THE CONTACT 
NONLINEARITY IN THE VICINITY OF 
VORTEX LINES 



The dynamics of vortex lines in Bose-Einstein conden- 
sates in a symmetric trap is governed by the GP equation 



indt^{v,t) = - A + — - 
' 2m 2 



^(r,t) 



(33) 



where a is the scattering length. Although we cannot 
solve this equation analytically, we can estimate the in- 
fiuence of the nonlinear term on the motion of vortex 
lines. Since the wave function vanishes on the vortex 
lines, the influence of the nonlincarity is weakening when 
one approaches the vortex line. On the other hand the 
interesting topological properties of the vortex lines mo- 
tion (reconnection, vortex creation, etc.) occur at small 
distances. Therefore, we restrict our analysis to distances 
much smaller (at least ten times) than the extension of 
the condensate. As an example we consider the inter- 
action of two perpendicular vortex lines separated by a 
distance d. The initial wave function is assumed in the 
form 

V'(r,t = 0) (34) 
= A{x + i{z - d/2)){y + i{z + d/2)) exp(-rVi^), 



where A = ViV7r-3/4i-7/2(3/2 + d7l6L4)-i/2 ig the 
normalization constant, L is the linear dimension of the 
condensate, and N is the number of atoms in the conden- 
sate. The ratio of the nonlinear part of the GP equation 
to the kinetic energy part for this wave function is 



32Ara e-'-'/^'(a;2 + (z - d/2)2)(y2 + + d/2)2)(a; + i{z - d/2)){y + i{z + d/2)) 



V^L3(3/2 + d4/i6L4)(8x4 + 28L^{x + iz)(y + iz) + d^{3L^ - r^) - 4r^{x + iz)(y + iz) + 2id{x - y)i5L^ 



(35) 
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At the center this ratio is equal to 



IV. CONCLUSIONS 



-All; 



Nacf 



r=0 



(36) 



Therefore it is smah, of the order of ^ = N{a/L){d/L)^. 
Taking the Unear dimension of the trap L = 5 x lO^'^ni, 
the vortex separation d = L/10, the number of atoms 
N = 10^, and the scattering length a — 5 x 10~^m, one 
obtains ^ = 10^'^. The same estimate of the ratio ( |35| ) 
is valid for all points lying at distances of the order of d 
from the center. Next, we have to estimate the role of the 
trap potential. The potential term modifies the shape of 
the wave function for times comparable with the period 
of the trap T w lO^^g, xhe characteristic time scale for 
the motion of vortex lines, as discussed in detail in is 
Tq = md^ /Ti. For the sodium atoms the value of this pa- 
rameter is To 3 X lO^'^s. Therefore, the role of the trap 
potential is negligible for small vortex separations and 
for times of the order of Tg. In Ref. Q we used the lin- 
ear Schrodinger equation to study the evolution of wave 
function that was initially given by Eq. (|3^). Vortex lines 
that initially were separated by a distance d crossed and 
reconnected after the time of the order of Tq. In view of 
our present analysis the same behavior should be found 
also in the evolution governed by the GP equation. The 
same arguments can be used in the case of more com- 
plicated vortex structures provided we restrict ourselves 
to distances small as compared to the linear dimension 
of the condensate and to times small as compared to the 
trap period. 



The aim of this study was to find out to what extent 
the topological properties of vortex lines motion carry 
over from linear to nonlinear Schrodinger equations. We 
analyzed two extreme cases of the nonlinear term that 
describes the mutual interaction of particles: the har- 
monic interaction (extremely long range) and the con- 
tact interaction (zero range). In the first case we were 
able to give a mathematical proof that the only change in 
the motion of vortex lines is an overall, time dependent 
displacement of the whole vortex line structure. In the 
second case we have shown that the features of the vortex 
lines motion that occur at small distances and for short 
times are governed mainly by the kinetic energy term in 
the Schrodinger equation. The influence of the nonlinear 
contact term under those conditions is negligible. Thus, 
many features of the vortex lines motion (especially those 
that involve a close approach of two vortex lines) are to 
a large extent universal. 
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